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Abstract. We compute the Parisi overlap distribution for paperfolding sequences. It 
turns out to be discrete, and to live on the dyadic rationals. Hence it is a pure point 
measure whose support is the full interval [—1, +1]. The space of paperfolding sequences 
has an ultrametric structure. Our example provides an illustration of some properties 
which were suggested to occur for pure states in spin glass models. 

1 Introduction 

In [TTJ the study of the Parisi overlap distribution for various classes of non-periodically 
ordered sequences was undertaken. 

The considered sequences were members of { — , +} z and their orbit closure forms typically 
a uniquely ergodic system, with a unique shift-invariant measure {i. This measure then 
can be a ground state for a translation-invariant interaction which one can construct 
[31 |2U [29], and the individual sequences will be the pure (extremal) ground states. 
Their overlap distribution gives the behaviour under the product measure of with itself, 
describing a two-replica system, written as [i x fjf of the overlap between two randomly 
(from this product measure) chosen, bi-infinite sequences a and a': 




(1) 



i=l,...N 
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Note that the overlap q aa > between two sequences is directly related to their Hamming 
distance dnip-, c') = 1 ~ q {"' ■ 

If we take two random sequences, each chosen according to the same shift-invariant 
measure, with probability one the above limit will exist, as follows from the ergodic 
theorem. The product measure of two ergodic measures, however, although shift-invariant, 
in general will not be ergodic, thus one may obtain different values for the overlap with 
positive probability. 

A simple example illustrating this is the symmetric measure which give equal weights 
to the two alternating sequences: 

fi a it(cr) = ^(5(+-... 1 a) + 5(- + ... 1 a)) (2) 

This measure is ergodic under the shift transformation, but the product measure of 
Halt with itself is not. For the overlaps it holds that with probability | the same sequence 
is chosen (overlap one) and with probability ^ two different sequences are chosen (overlap 
minus one). Thus in this simple example 

p(g) = 1) + % -1)) (3) 

In more general situations the overlap can take any value in the interval [—1,-1-1], 
and the induced measure, which is the overlap distribution, describable by a probability 
distribution function on this interval, can be discrete ( that is, it is a pure point measure) 
or nondiscrete, and thus have continuous components. 

In particular, it was found in [11] that whenever the atomic ("diffraction") spectrum of 
the unique translation-invariant measure on the sequences does not contain a pure point 
component, the overlap distribution is trivial; this applies in particular to all weakly 
mixing systems (in which case their product measures are ergodic), as well as to the 
(Prouhet-(Thue-))Morse system. 

On the other hand for the Fibonacci sequences the overlap distribution was found to 
contain a(n absolutely) continuous part. In fact, in explicit form it is [UJ 

p(q)dq = (2 7 - 1)%, 1 - 4(1 - 7 )) + l -l {l ^ llx] {q)dq (4) 

where 7 = : 2 7 = . 

' l+Vo 

We notice here that the same argument which was used in [TT] applies to more general 
Sturmian sequences, which thus also have a continuous part in their overlap distribution. 
The reason is that these are also rotation sequences, see e.g. [IT]. They also have recently 
been studied as ground states of some explicit interaction functions on lattice systems in 

&)■ 

Moreover, it was found that the Toeplitz (or period-doubling) sequences give rise to 
an overlap distribution which is concentrated on a countable number of values, with 1 
as its only limit point. Also, the set of these Toeplitz sequences has a tree (ultrametric) 
structure, and the overlap distribution could be obtained via this structure. 

Here we show that the paperfolding sequences (see [5J [H] and references mentioned 
there), which also display such an ultrametric structure, give rise to an overlap distribution 
living on a dense countable number of points. As such the paperfolding system, which, in 



3 



contrast to the Toeplitz system, is symmetric under the spin-flip (plus-minus) symmetry, 
is even closer to what is expected in Parisi's replica symmetry breaking theory of the 
Sherrington-Kirkpatrick model [18j|27]. I n Parisi's theory for the SK model, this countable 
overlap distribution is supposed to show up for a fixed disorder, while averaging over the 
disorder distribution is supposed to result in a continuous piece in the overlap distribution. 
Often the countable overlap distribution is interpreted as being related to or implying a 
countable number of pure states; however there seems not to be much justification for 
this 0. 

Conceptually, our example illustrates a number of points: 

Newman and Stein [2U [221 E3J [2U 1251 ED] have shown that self-averaging arguments, 
based on the use of the ergodic theorem, apply to many finite-dimensional, disordered 
(spin-glass) models (although not to the equivalent-neighbour Sherrington-Kirkpatrick 
model). In particular they find that for systems with a spatial structure, in contrast 
to the Parisi theory's predictions for the Sherrington-Kirkpatrick model [HI [27], overlap 
quantities cannot depend on the realization of the disorder parameter, that is, they are 
"self-averaging". For example, for Edwards- Anderson-type spin-glass models there is no 
difference between a typical (Parisi SK-prediction: supported on countably many values, 
dense) and an averaged (Parisi SK-prediction: continuous) overlap distribution. 

As mentioned before, for many nonperiodic systems, the individual (nonperiodic) se- 
quences are ground states for some translation invariant deterministic interaction. Thus 
they are often used as models for non-periodic (quasi-) crystals, [21 [2E1 1221 EH EH HEH ED 
[33] . It therefore becomes natural to consider various questions naturally arising in statisti- 
cal mechanics about such systems, whether about their spectral properties [H [E2J EES HE] 
or, as was done in [TH |10] and also here, about their overlap distributions. Typically, 
there are uncountably many such ground states, and as mentioned before, they can have 
nontrivial overlap distributions; somewhat surprisingly even without the presence of any 
disorder. 

It turns out that: 

• There is no general connection between the number of pure states -uncountable- and 
the number of possible overlap values -which can be one, countably or uncountably 
many-. 

• Triviality, non-triviality, and ultrametricity of the Parisi overlap distribution all can 
be realized without any role for the disorder. 

• A dense discrete overlap distribution, with a hierarchical structure, can be realised 
(by the paperfolding system). 

2 Main Result 

For background on the paperfolding system we refer to [U [2] , and references mentioned 
there. 

A paperfolding sequence can be constructed as follows: 

In step 1, we choose k\ to be a number from the set {1, 2, 3, 4}. 
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Now we fill all sites of the form k\ + 4Z with pluses, and all sites of the form k\ + 2 + 4Z 
with minuses. 

Thus we have the period-4 structure covering half the sites: 
+ .- . + .- . + .- . + .- .... 

In step 2, we again choose &2 to be a number from the set {1, 2, 3, 4}. 

Now we occupy all sites of the form fc 1 + l+2A;2+8Z with pluses and all sites k 1 +5+2k2+8% 

with minuses. 

We repeat this step and in the mth step we again occupy one quarter of the remaining 
empty sites with pluses, periodically with period 2 x 2 m , and the quarter of the sites at 
distance 2 m from the pluses we occupy with minuses, periodically with the same period. 
Thus in each step one half of the remaining sites are filled. 

After making a choice for k m , for all m G N, all sites will be filled, and we have 
a paperfolding sequence. Notice that in every step the uncountable set of sequences is 
divided into four times as many subsets as in the previous step, which generates the 
hierarchical (ultrametric) structure of the uncountable set of sequences. 

With probability | we find the overlap value q = 0, which happens if in step 1 one 
"odd" and one "even" sequence was chosen, that is, k\ — k[ is odd. With probability | one 
finds overlap values q equal to plus or minus = respectively, which happens if in the first 
step either the same (ki = k[), or the opposite (\ki — k[\ =2) choice was made, and in step 
2 a different even-odd choice for &2 and k' 2 . We repeat this argument, and we find that 

with probability 22 l +1 one finds overlap values ^ 2p 2 t , with p = — 2 n_1 , , 2™ _1 — 1 [T4] . 

Thus in the end, the overlap distribution becomes concentrated on the dyadic rationals. 
Summarizing we have proven: 

Theorem 1 The overlap distribution of paperfolding sequences is given by: 

OO ' 1 / 1 \ " 

Hi) ^) (5) 

n=0 m v / 

where the integers n, m must satisfy the following conditions: 

• if n = then m = 0, 

• if n > then m is odd and \m\ < 2 n . 
Comments: 

Properties which are common between the period-doubling and the paperfolding sys- 
tem are: 

1) Limit periodicity. The sequences are a countable union of periodic subsequences of 
pluses and minuses. This implies a countable number of possible values for the overlap. 

2) There is a tree (hierarchical, ultrametric) structure in the set of sequences- the pure 
states-. As the overlap between two sequences depends on the level of the tree at which 
a different even-odd choice is made, the overlap distribution reflects this. 

3) Both paperfolding and period-doubling sequences form model sets, obtained by a 
cut and project scheme, with an internal space of 2-adic numbers. [61 El [32] . This agrees 
with the existence of a connection between Parisi's replica-symmetry-breaking ideas and 
p-adic numbers as was suggested in [2^] . 
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However, in contrast to the period-doubling case, the set of paperfolding sequences is 
spin-flip symmetric, and the set on which the overlap distribution is concentrated now 
lies dense in the interval [— 1, +1]. 

Acknowledgements: A. CD. v.E became interested in these issues during an earlier collab- 
oration with Bert Hof and Jacek Mi§kisz. He also acknowledges very useful discussions 
with Michael Baake, Jacek Miejdsz and Dimitri Petritis. We also thank Michael Baake 
for some helpful advice on the manuscript. 



References 

[1] V. Anagnostopoulou, K. Diaz-Ordaz, O. Jenkinson and C. Richard, newblock Stur- 
mian maximizing measures for the piecewise-linear cosine family. Dresden preprint, 
2010. 

[2] J. -P. Allouche and J. Shallit Automatic sequences. Cambridge University Press, 
Cambridge, 2003. 

[3] S. Aubry. Weakly periodic structures and example. J.Phys. (Paris) Coll., C3-50:97- 
106, 1989. 

[4] M. Baake and M. Hoffe. Diffraction of random tilings:Some rigorous results. J. Stat. 
Phys. 99:219-261, 2000. 

[5] M. Baake, R. V. Moody, C. Richard, B Sing. Which distributions of matter diffract ?- 
Some answers. Quasicrystals: Structure and Physical Properties: Ed. H.-R. Trebin, 
Wiley- VCH p.188-207, arXiv:math-ph/0301019, 2003. 

[6] M. Baake, R. V. Moody and M. Schlottmann. Limit-(quasi)periodic point sets as 
quasicrystals with p-adic internal spaces. J. Phys. A, Math. Gen. 31:5755-5765, 1998. 

[7] M. Baake, R. V. Moody. Weighted Dirac combs with pure point diffraction. J. Reine 
und Angewandte Mathematik 573, p61-94, arXiv:math.MG/0203030, 2004. 

[8] E. Bolthausen. In Spin glasses, Eds.E. Bolthausen and A. Bovier. Springer LNM 
1900. 2007. See in particular pl6. 

[9] M. Dekking, M. Mendes-France, A. van der Poorten FOLDS! (1,11 and III) The 
Mathematical Intelligencer, 4, p. 130-138, p. 173-181, p. 190-195, 1982. 

[10] A. C. D. van Enter. On the set of pure states for some systems with non-periodic 
long-range order. Physica A 232:600-607, 1996. 

[11] A. C. D. van Enter, A. Hof and J. Miejdsz. Overlap distributions for deterministic 
systems with many pure states. J. of Phys. A, Math. Gen. 25:L1133-1137, 1992. 

[12] A. C. D. van Enter and J. Miejdsz. How should one define a (weak) crystal? J. Stat. 
Phys. 66:1147-1153, 1992. 



6 



[13] C. Gardner, J. Miejdsz, C. Radin and A. C. D. van Enter. Fractal symmetry in an 
Ising model. J. Phys. A, Math. Gen., 22:L1019-1023, 1989. 

[14] E. de Groote. The overlap distribution of paperfolding sequences. Groningen bachelor 
thesis, 2010. 

[15] J.-Y. Lee, and R. V. Moody. Lattice substitution systems and model sets. Discrete 
and Computational Geometry 25:173-201, 2001. 

[16] J.-Y. Lee, R. V. Moody and B. Solomyak. Pure point and dynamical diffraction 
spectra. Ann. Inst. Poincare, Theor. Math. Phys. 3:1003-1018, 2002. 

[17] K. Lii and J. Wang Construction of Sturmian sequences. J. Phys. A: Math. Gen. 
38:2891-2897, 2005. 

[18] M. Mezard, G. Parisi and M. A. Virasoro. Spin Glass Theory and Beyond. World 
Scientific, Singapore, 1987. 

[19] J. Mi§kisz. Quasicrystals. Microscopic models of non-periodic structures Leuven 
Lecture Notes in Mathematical and Theoretical Physics, Vol. 5, 1993. 

[20] C. M. Newman. Topics in Disordered Systems. ETH Lectures in mathematics, 
Birkhauser, Basel, Boston, Berlin 1997. 

[21] C. M. Newman and D. L. Stein. The metastate approach to thermodynamic chaos. 
Phys. Rev. E 55: 594-5211, 1997. 

[22] C. M. Newman and D. L. Stein. Thermodynamic chaos and the structure of short- 
range spin glasses. Mathematical Aspects of Spin Glasses and Neural Networks 243- 
287, eds A. Bovier and P. Picco, Birkhauser, Boston-Basel-Berlin, 1998. 

[23] C. M. Newman and D. L. Stein. Ordering and Broken Symmetry in Short-ranged 
Spin Glasses. arXiv:cond-mat/0301403 J. Phys., Cond. Matt, 15, R1319-R1364, 
Topical Review, 2003. 

[24] C. M. Newman and D. L. Stein. The state(s) of replica symmetry breaking: Mean 
field theories versus short-ranged spin glasses. (Formerly known as "Replica Symme- 
try Breaking's New Clothes") J. Stat.Phys. 106: 213-244, 2002. 

[25] C. M. Newman and D. L. Stein. Distribution of Pure States in Short-Range Spin 
Glasses Int. J. Mod. Phys., B24, 2091-2106, 2010. 

[26] G. Parisi and N. Sourlas. P-adic numbers and replica symmetry breaking. Eur. J. 
Phys., B14, 235-242, 2000. 

[27] G. Parisi and M. Talagrand. On the distribution of the overlaps at given disorder. 
C.R Acad. Sci. Paris, Ser I, 339, 303-3006, 2004. 

[28] C. Radin. Disordered ground states of classical lattice models. Rev. Math. Phys., 
3:125-135, 1991. 



7 



[29] C. Radin. Low temperature and the origin of crystalline symmetry. Int. J. Mod. 
Phys. B, 1:1157-1191, 1987. 

[30] D. Ruelle. Do turbulent crystals exist? Physica, 113A:619-623, 1982. 

[31] D. Schechtman, I. Blech, D. Gratias and J. W. Calm. Metallic phase with long-range 
orientational order and no translation symmetry. Phys. Rev. Lett., 53:1951-1953, 
1984. 

[32] M. Schlottmann. Generalized model sets and dynamical systems. : In Directions in 
Quasicrystals, 143-159, eds M. Baake and R. V. Moody, CRM Monograph Series, 
vol 13, AMS, providence, RI, 2000. 

[33] M.Senechal. Quasicrystals and geometry. Cambridge University Press, 1995. 



8 



